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We present results of searches for B-meson decays to K + it~-k + with the BABAR detector. With 
a data sample of 61.6 million BB pairs, we measure the branching fractions and 90% confidence- 
level upper limits averaged over charge-conjugate states (the first error is statistical and the sec- 
ond is systematic): B(B + -> K *°(892)7r+) = (15.5 ± 1.8^) x 10" 6 , B(B+ f (980)K+J -> 
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TT+iY-) = (9.2 ± 1.2±|;J) x 10~ 6 , B(B+ -> D°ir + y D° -» AT+tt-) = (184.6 ± 3.2 ± 9.7) x KT 6 , 
B(B + -> p°(770)K + ) < 6.2 x 10" 6 and B{B + -> A" + 7r-7r + non-resonant) < 17 x 10" 6 . 



PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 

The study of charmless hadronic B decays can 
make important contributions to our understanding of 
hadronic decays and CP violation in the Standard 
Model 0. Branching fraction predictions for B meson 
decays to Pseudoscalar- Vector final states have recently 
been calculated using QCD Factorisation and SU(3) fla- 
vor symmetry models L3- Li| • The measurement of B + 
meson [j| decays to the final state K + -k~-k + via interme- 
diate resonances can be used to search for weak phases 
and direct CP violation. The signal can be used to ex- 
amine the light-meson mass spectrum 0, . The charm 
state Xco might be sensitive to the angle 7 of the Unitarity 
Triangle through interference with the non-resonant com- 
ponent producing an observable charge asymmetry ||; 
the branching fraction also constrains some models for 
charmonium hybrid production |9j]. 

The data used in this analysis were collected at the 
PEP-II asymmetric e + e~ storage ring with the BABAR 
detector [10| . The BABAR detector consists of a five-layer 
silicon tracker, a drift chamber, a new type of Cherenkov 
detector [ill ] , an electromagnetic calorimeter and a mag- 
net with instrumented flux return. The data sample has 
an integrated luminosity of 56.4 fb _1 collected at the 
T(4S) resonance, which corresponds to (61.6 ±0.7) x 10 6 
BB pairs (N B g). We assume that the T(4S) decays 
equally to neutral and charged B meson pairs. 

The K + Tt~Tr + phase space can be represented in a 
Dalitz plot, in which the many resonant B decay modes 
form overlapping bands and interference occurs where the 
bands overlap 12]. As a consequence, the whole Dalitz 
plot should be considered before assigning a branching 
fraction to a specific mode. However, the data sample 
is not large enough for a full Dalitz plot fit to be ef- 
fective. In this analysis, the Dalitz plot is divided into 
eight regions, reflecting as much as possible the known 
decay modes. We first determine the yields in these re- 
gions using a maximum- likelihood fit. We then interpret 
the yields as branching fractions, assuming a particular 
collection of quasi-two-body decay modes. We evaluate 
the systematic uncertainty due to the particular choice 
of decay modes and the effect of interference between the 
different contributions. 

The regions are defined in ttik-k and m W7r , the in- 
variant masses of the neutral Kir and 7T7t systems, as 
given m Table [Q and illustrated in Fig. [T] Re gion I is 
expected to be dominated by iC* (892)7r + . Region II 
could have contributions from several higher K* reso- 
nances. Region III is dominated by the production of 
D°tt + . The high branching fraction for this mode al- 
lows it to be used to correct for differences between data 
and simulated events and to evaluate systematic uncer- 




FIG. 1: A Dalitz plot showing the defined regions. 



tainties. Regions IV and V are expected to be domi- 
nated by p°(770)K+ and f (980)K+, respectively. The 
resonance contributions to Region VI are not known a 
priori. The area where these regions intersect the D° 
band, 1.8 < tuk-k < 1.9GeV/c 2 , is excluded. Region 
VII could contain higher mass charmless and charmo- 
nium resonances, as well as a non-resonant contribution 
that extends across the whole Dalitz plot. Region VIII is 
dominated by XcqK + . This channel is vetoed from other 
regions using 3.355 < < 3.475 GeV/c 2 . 

TABLE I: Regions of the K + it~-k + Dalitz plot and signal 
yields obtained (the first error is statistical and the second 
is systematic). The D° band, 1.8 < ttikti < 1.9GeV/c 2 , is 
excluded from all regions except region III, and the Xco band, 
3.355 < rriTTTT < 3.475 GeV/c 2 , is excluded from all regions 
except region VIII. 



Region 


Selection Criteria (GeV/c 2 ) 


Signal Yield 


I 


0.816 < m K n < 0.976 , > 1.5 


161 ± 18 ±4 


II 


0.976 < mxvr < 1-8 , m n7: > 1.5 


405 ± 28 ± 13 


III 


1.835 < m K Tt < 1-895 


3755 ±66 ±11 


IV 


0.6 < m W7T < 0.9 


66 ± 15^ 


V 


0.9 < m vv < 1.1 


179 ± 19 ± 5 


VI 


1.1 < m wlT < 1.5 


126 ± 19 ± 5 


VII 


rriKir > 1-9 , > 1.5 


133 ± 23122 


VIII 


rriK-n > 1-9 , 3.37 < m^^ < 3.46 


26 ± 6 ± 1 



Candidate B mesons are formed by combining three 
charged tracks, where each track is required to have at 
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least 12 hits in the drift chamber, to have transverse 
momentum of at least 100 MeV/c and to be consistent 
with originating from the beam-spot. Charged pions and 
kaons are identified using energy loss (dE / dx) measured 
in the tracking system and the Cherenkov angle and num- 
ber of photons measured by the Cherenkov detector. Pi- 
ons are required to fail the kaon selection. The efficiency 
of selecting kaons is approximately 80%, while the proba- 
bility of misidentifying pions as kaons is below 5%, up to 
a laboratory momentum of 4.0 GeV/c. Pions are also re- 
quired to fail an electron selector which uses dE/dx, the 
energy to momentum ratio and the shape of the calorime- 
ter signal. Over 99% of pions from the signal decay pass 
this requirement. 

Signal decays are identified using two kinematic vari- 
ables: AE, the difference between the center-of-mass 
(CM) energy of the B candidate and \fsj2, where s/s is 
the total CM energy; and the beam-energy-substituted 
mass toes = \J (s/2 + p { ■ p B ) 2 /Ef - p|, where p B is 
the momentum of the reconstructed B candidate and 
(Ei,pi) is the four-momentum of the initial e + e~ sys- 
tem. The AE and toes distributions for signal events 
have widths of 20 MeV and 2.7 MeV/c 2 , respectively. We 
require 5.22 < to es < 5.29GeV/c 2 and \AE\ < 0.1 GeV 
for events entering the fit. Events with 0.1 < \AE\ < 
0.3 GeV are used for continuum background characteri- 
sation as described below. 

A very small proportion of events, fewer than 4%, have 
two or more candidates that pass the above requirements. 
For these events a single candidate is selected at random, 
so as not to bias the fit distributions. This random se- 
lection has a minimal impact on the efficiency, and any 
systematic uncertainty due to this effect is negligible. A 
single candidate per event is similarly selected for the 
data with 0.1 < \ AE\ < 0.3 GeV used in continuum back- 
ground characterisation. 

Continuum light-quark and charm production is the 
dominant source of background. This is suppressed us- 
ing two event-shape variables. The first is the cosine 
of the angle 0t between the thrust axis of the selected 
B candidate and the thrust axis of the rest of the event. 
For continuum background, the distribution of | cos 6t\ is 
strongly peaked towards unity whereas the distribution 
is uniform for signal events. We require |cos#t| < 0.9. 
The second event-shape variable is a Fisher discriminant 
(F) ^3 ■ F° r -F we use a linear combination of the cosine 
of the angle between the B-candidate momentum and 
the beam axis, the cosine of the angle between the B- 
candidate thrust axis and the beam axis, and the energy 
flow of the rest of the event into each of nine contigu- 
ous, concentric, 10° cones around the thrust axis of the 
reconstructed B [lij . 

There are also B-decay backgrounds, mainly four-body 
decays, and three-body decays with one or more parti- 
cles misidentified. These backgrounds are studied us- 
ing Monte Carlo simulations (MC). They are reduced 



by the particle identification selections and by exclud- 
ing events containing J/i/j or ip(2S) decays to l + l~ with 
vetoes 2.97 < m„ < 3.17 GeV/ c 2 and 3.56 < % < 
3.76 GeV/c 2 , however some -B-decay backgrounds remain 
and must be accounted for. For backgrounds contribut- 
ing only a few events to the maximum-likelihood (ML) 
fit, the estimated contribution is subtracted from the fi- 
nal signal yield with a systematic uncertainty to account 
for the unknown probability of the background to be se- 
lected as signal in the fit. Larger backgrounds are pa- 
rameterized in the ML fit as described below. These 
were B + — > D°7r + , D° — > K + ir~ir° in regions II and VII, 
B+ -► D°p+ (770) with D° -> K+ir- and p+ -> tt+tt 
in region VII and B + — > rj'K + with rf — > p° (770)7, 
p° — > 7r + 7r~ in regions IV and V. 

We form probability density functions (PDFs) with pa- 
rameters a for the three variables (x) toes, AE, and T in 
each region. We find the correlations among these vari- 
ables to be negligible; accordingly, for each hypothesis I 
(signal, continuum background, and B background), we 
form a product Pi = Pi.ttissPi^ePi^ that models that 
hypothesis. The likelihood for an event j is the sum over 
the M hypotheses of the products Pi , with each product 
weighted by the number of events (to be determined), m, 
for that hypothesis. A product over the N events in the 
data sample of the event likelihoods along with a Poisson 
factor forms the likelihood function: 

(M \ N / M \ 
- e ni J n (e ni Pi & 3) j ■ w 

This likelihood is maximized to obtain ni for the sig- 
nal and continuum background components; n; for the 
B-background component is fixed to an estimate of the 
contribution from MC. The parameters of the signal and 
J5-background PDFs are determined from MC and held 
fixed in the final fit. The continuum background param- 
eters for T are fixed from the data with 0.1 < \AE\ < 
0.3 GeV. The continuum background parameters for AE 
and Toes are left free in the final fit. Each toes PDF is a 
Gaussian distribution for signal, and an ARGUS thresh- 
old function [l^ for continuum background. Each AE 
PDF is a sum of two Gaussian distributions with equal 
means for the signal and a first-degree polynomial for 
the continuum background. The signal and background 
T PDFs are Gaussian distributions with different widths 
above and below the mean. For the B-background pa- 
rameterizations, signal or continuum shapes are used de- 
pending on the nature of the background. The projec- 
tions of the Region I data in m^s, AE and T and the re- 
sults of the fit are shown in Fig. El demonstrating a clear 
signal. The data in the plots pass a selection on the per- 
event signal-to-background likelihood ratio formed from 
the other two fit variables, which has been optimized to 
give the greatest significance to the signal. 

The signal yields for the regions of the Dalitz plot 
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Fisher Discriminant 

FIG. 2: Projection plots in rriES, Ai? and T , for the data in 
Region I. The superimposed curve is a projection of the full 
fit with the background component shown as a dotted line 
and, for T, the signal component shown as a dashed line. 



TABLE II: A summary of the model used to calculate branch- 
ing fractions. Alternative lineshapes are given in parentheses. 
The masses and widths are taken from the Review of Particle 
Physics 0. 



Resonance 


Lineshape 


Mass 


Width 






(MeV/c 2 ) 


(MeV/c 2 ) 


iT u (892) 


BW 


896.10 ±0.27 


50.7 ±0.6 


is:o* (1430) 


BW (LASS[1£]) 


1412 ± 6 


294 ± 23 


D° 


BW 


1864.5 ± 0.5 





p°(770) 


Blatt-Weisskopf[17j 


769.0 ± 0.9 


150.9 ±1.7 


/o(980) 


BW (Flatte[18j) 


980 ± 10 


70 ±30 


/ 2 (1270) 


BW 


1275 ± 12 


185 ± 30 




BW 


3415.1 ±0.8 


16.2 ± 3.2 


non-resonant 


flat 


all masses 





are shown in Table [fl The systematic uncertainty 
arises from the PDF parameters and from B-background 
subtraction. We find that all yields have a signif- 
icance greater than five standard deviations, where 
J2ln(C m ax/ '£ m „(„ sijiii]F o)) is used as an estimator of the 
significance. 

We calculate the branching fractions from the mea- 
sured yields, taking into account the overlapping nature 
of the resonances, using B = MT 1 Y/N BB where Y is a 
vector of the yields in each Dalitz region and B is a vector 
of branching fractions. M is the efficiency matrix where 
My is the probability that an event arising from the con- 
tribution dominating region i will be found in region j. 
The elements of M are estimated using MC including 
small corrections for differences in tracking and in parti- 
cle identification efficiencies between MC and data, and 
differences between MC and our resonance model. 

In our resonance model we assume one dominant con- 
tribution per region. The contributions for the chosen 
model are given in Table [H] For Regions II and VI 
where the main contributions are not known a priori, we 



TABLE III: The measured branching fractions and uncertain- 
ties. The first uncertainty is statistical, the second includes 
the systematic uncertainties from the yields and the efficien- 
cies, the third is the model uncertainty, and the fourth is the 
uncertainty due to interference. 



Channel 
jr°(892)7r + 

"higher K*°"ir+, K*° -> K+ty- 
D° 7T+, D° -» K+7V- 
p°(77Q)K+ 

"higher f"K+, f -> tt+tt^ 
Non-resonant 

Xc0K + , XcO -> 7T + 7T~ 



BF xlO' 6 

15.5 ±1.8 ±1.1^;| ±0.9 
25.1 ± 2.0 ± 2.9±s * ± 4.9 
184.6 ±3.2 ±9.7 
3.9±1.2+0; 3 6 +° 3 ; 3 2 ±1.2 

9.2±1.2±0.6±i;|±1.6 
3.2 ± 1.2 ± 0.5±1;1 ± 1.5 
5.2 ± 1.9±S;l±r.l ± 6.4 
1.5 ±0.4 ±0.1 



take the dominant contributions to be ifg°(1430)7r + and 
/2(1270)if + respectively. However, we quote branch- 
ing fractions for B + —> "higher K* 0, 'ir + where "higher 
K*°" means any combination of Kf (1430), Kf (1430) 
and *rf(1680) and B+ ->■ "higher f"K+, where "higher 
/" means any combination of / 2 (1270), /n(1370) and 
/2(1430). Non-relativistic Breit-Wigner line shapes are 
used for all channels except for the broad p(770) res- 
onance, where we use a relativistic Breit-Wigner line 
shape with Blatt-Weisskopf damping |lj| ■ 

We evaluate systematic uncertainties on the branching 
fractions taking into consideration uncertainties on res- 
onance parameters and alternative line shapes given in 
Table llll We also include the effect on all branching frac- 
tions that in Regions II and VI, the dominant contribu- 
tion could be any combination of a number of resonances. 
Also included in the systematic uncertainties is the pos- 
sibility that the yield measured in Region VII is from a 
resonant component and does not extend into the other 
regions. Uncertainties on the branching fractions due to 
the interference between the resonances are evaluated by 
generating many simulated B + — > K + t:^t: + Dalitz plots 
with all the contributions having random phases and ob- 
serving how the interference between the contributions 
affects the measured branching fractions. The branching 
fractions and uncertainties of intermediate resonances are 
given in Table ITTTl 

Fig. |3 shows the Dalitz plot for data events within 
a signal region, 5.2715 < ttt-es < 5.2865 GeV/c 2 , that 
have a the per-event signal-to-background likelihood ra- 
tio formed from the AE and J- PDFs, greater than 5. 
Both signal and background events appear in the plot. 
The D° 7r + signal is the narrow band in the mjf, spec- 
trum. To illustrate the expected background distribu- 
tion, events passing the same likelihood selection but hav- 
ing a value of jties between 5.25 < ities < 5.26GeV/c 2 
are also shown. The size of this sideband is chosen to con- 
tain approximately the expected number of background 



7 



IS* 



10 



(3 

-20 • 



15 20 25 
m 2 K .JGeVlc 2 ) 2 



10 15 20 25 

mj^XGeV/c 2 ) 2 



FIG. 3: Dalitz plots showing (left) the observed distribution 
in a signal jties region (defined in the text) and (right) the dis- 
tribution for continuum background from the tties sideband. 
Vetoes reject events with J/xj) and ip{2S). 



events that will enter the signal-region plot. 
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FIG. 4: Background-subtracted and efficiency-corrected pro- 
jections of the Dalitz plot in mx, and m,™- 

Fig. 0] shows background-subtracted, efficiency- 
corrected projections of the two-body invariant mass 
spectra, tuk-w and m^n, from 0.6GeV/c 2 to 1.8GeV/c 2 
and from 0.2GeV/c 2 to I.5GeV/c 2 respectively. The sig- 
nal events and background distributions are obtained by 
the same method as for Fig. |3| and again the D°, J/^ 
and ^(2S) vetoes are applied. Peaks at the if* (892) 
and /o(980) masses are clearly visible. 

The invariant mass m^i or m V7r , and helicity angle be- 
tween the resonance decay and flight directions, Oh, are 
not used in the likelihood fit. However, to illustrate our 
findings, we show, in Fig. resonant mass and cos Oh 
projections for Regions I, IV and II after background 
subtraction and efficiency corrections. Figs.E{a-d) have 
been overlaid with the distribution of the expected dom- 
inant resonance: Breit-Wigner line shapes for the mass 
distributions in Figs. |3fa) and OJc), cos 2 H for the 
if* (892) cos 9h distribution in Fig. GJb), and a uniform 
distribution for the scalar / (980) cos Oh distribution in 
Fig-EJd). There is good agreement between the overlaid 
and observed distributions indicating that the expected 
resonances are indeed dominant in these regions. The 
Jo (980) cos Oh distribution suggests a linear dependence 
that is most likely due to interference with the vector 
p(770), which is taken into account in our interference 
systematic uncertainty. 




0.85 0.9 0.95 

m Kii (GeV/c 2 ) 



-1 -0.5 0.5 1 
Cos 9 H 




^200 



rnoo- 



~ 





-u+' 





1.2 1.4 1.6 1.8 
m K ,(GeV/c 2 ) 



CM 

O200 


-\ 





\ 




W 

100 

LU 




+ - 




+ + + 

















1 


-0.5 





0.5 



Cos 9 U 



FIG. 5 

cos 8h 



Projection plots of the two-body invariant mass and 
for, from top to bottom, Regions I, V and II. 



We can see in Fig0J there is a large signal in the region 
1.1 < niK-K < 1.4GeV/c 2 (Region II). This is shown in 
more detail in Fig.[S{e), resonant mass, and Fig.OJf), the 
cos Oh projections for Region II. The complex behavior 
of this signal is similar to that observed by LASS 0] 
and precludes an interpretation as a single resonance. 

In conclusion, we have made branching fraction mea- 
surements, summarized in Table IIIII for a number of 
charm and charmless B decay channels with the final 
state K + tt~it + . This analysis has taken into account 
the uncertainty in the knowledge of the nature and pa- 
rameterization of the intermediate resonances on all the 
branching fractions assuming a non-resonant contribu- 
tion with kinematics defined by phase space. The re- 
sults also take account of the unknown levels of in- 
terference between the different contributions. The 
B + — ► D°tt + and B + -±_x c oK + results agree with previ- 
ous measurements [HEf. The B+ K*°(892)ir+ H3 
and B + — > /o(980)A + branching fractions are consistent 
with, and more precise than, previous measurements • 
The B + — > A*°(892)7r+ result is significantly higher than 
predicted by many factorization models [3( . The observa- 
tion of the decay B + — > /o(980)A + provides hints about 
the nature of the / (980) Q. A large signal is seen for 
B + ^"higher K*°"tt + where "higher A* " means any 
combination of A,* o (1430), A 2 *°(1430) and A7°(1680). 

We also give 90% confidence-level upper limits for 
the branching fractions of the following channels: 
B(B+ p°(770)A+) < 6.2 x 10" 6 , B(B+ K+ir-% + 
non-resonant) < 17 x 10~ 6 , B(B+ -> "higher /"A+) < 
12 x 10~ 6 . The tight limit on the non-resonant compo- 
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nent means that its 7-dependent interference with the 
XcoK final state will be very hard to measure. 
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